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Abstract 



The dynamics of small perturbations around sphaleron and black hole solutions in 
the Einstein- Yang-Mills theory for the gauge group SU (2) is investigated. The perturba- 
tions can be split into the two independent sectors in accordance with their parity; each 
sector contains negative modes. The even-parity negative modes are shown to correspond 
to the negative second variations of the height of the potential energy barrier near the 
sphaleron. For the odd-parity sector, the existence of precisely n (the number of nodes 
of the sphaleron solution gauge field function) negative modes is rigorously proven. The 
same results hold for the Einstein- Yang-Mills black holes as well. 



^ Published in Geometry and Integrable Models, Eds. P.N. Pyatov & S.N.Solodukhin, World Scientific 
1996, pp. 55-77. 



1 Introduction 



Soon after the discovery of the Bartnik-McKinnon sohtons and black holes in the 
SU{2) Einstein- Yang-Mills (EYM) theory, it was shown that these interesting objects 
are unfortunately unstable 0. However, as has been realized later [Q, the property 
to be unstable is in fact fairly important in the case, since it allows for the sphaleron 
interpretation 0,0 of the regular BK objects. 

Usually, sphalerons are related with the existence of non-contractible loops in the 
configuration space of the theory 0, f^]- For each such a loop, one finds a maximal value 
of the potential energy and then minimizes the result over all loops; this gives the energy 
of a static saddle point field configuration. Notice that the configuration space has to 
be compact in this approach, otherwise the existence of non-contractible loops does not 
ensure the existence of saddle points of the energy The difference between the two 
cases is perfectly illustrated in the two figures in famous Manton's paper (see Fig.l and 
Fig. 2 in [Q). Although for the non-compact case the saddle points may also exist (see 
Fig. 2 in [^]), these points have more then one directions of instability, and thus can not 
be constructed in the same manner is in the compact case, such that their existence 
is not directly related to the existence of non-contractible loops. The use of loops in 
the non-compact case may have another meaning: for theories with vacuum periodicity, 
the non-contractible loop through the saddle point (if any exists) can be transformed 
into a path interpolating between distinct vacua. This shows that the saddle point field 
configuration relates to the top of the potential barrier between the vacua thus justifying 
the usage of the name "sphaleron" . 

It is not widely understood that the EYM sphalerons correspond to the "non- 
compact" case. In the EYM theory, the maximal value of energy for each non-contractible 
loop can be reduced to zero by a smooth deformation of the loop , such that the exis- 
tence of the loops does not ensure the existence of the solutions (otherwise the existence 
proof for the BK solutions - |10| would just follow the same line as that for the stan- 
dard sphaleron |TT[). Instead, in order to illustrate the sphaleron nature of the solutions, 
the loops can be used as paths connecting the solutions with the neighboring vacua 
(see below). Manifestly, the difference between the standard Yang-Mills-Higgs (YMH) 
sphaleron and the EYM sphalerons results in the different numbers of their unstable 
modes. The YMH sphaleron has one and only one negative mode in the odd-parity 
spherically-symmetric perturbation sector. In the even-parity sector, the solution has no 
negative modes, such that the energy of the solution indeed defines the minimal height 
of the barrier. The EYM sphalerons have negative modes both in the odd parity and 
in the even-parity sectors, which means that they do not specify the minimal possible 
height of the barrier (the latter is zero) |13|. 

The present paper concerns with the investigation of the dynamics of spherically- 
symmetric perturbation modes around EYM sphalerons and black holes. The perturba- 
tions decouple into the two independent groups in accordance with their parity. The even 
parity modes were studied previously in |]^ , it has been found numerically that the n-th 
soliton or black hole solution has precisely n negative modes of this type. Here we explic- 
itly show that the negative modes of this type are related to the decrease of the height of 
the potential barrier around the sphaleron |jl3[ . For the odd-parity modes, the existence 
of at least one such a mode for all known SU (2) EYM solitons and black holes has been 



established in M, WM, as well as for the "generic" solitons in the EYM theory for any 
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compact gauge group [ffSl . For the lowest n regular BK solutions, numerical analysis has 



revealed the existence of n negative modes of this type too [|l6l. Here we give a complete 
proof of the following statement: for all known regular and black hole non-Abelian EYM 
solutions there exist precisely n negative modes in the odd parity spherically-symmetric 



perturbation sector . For the even-parity modes the corresponding proof is unknown, 
however the results of the numerical analysis suggest that there are n negative modes in 
this case too, such that the total number of negative modes for the EYM solitons and 
black holes is 2n. 

The rest of the paper is organized as follows. In Sec. 2 the basic notations and 
equations are introduced. Sec. 3 contains an explicit construction of vacuum-to-vacuum 
paths through the EYM sphalerons. In Sec. 4 the basic perturbation equations are given 
for the both parity types. For the odd-parity perturbation modes, the transformation to a 
one-channel Schrodinger problem is given in Sec. 5, the regularity of the effective potential 
is shown in Sec. 6. A zero-energy bound state in the problem is explicitly presented in 
Sec. 7, the corresponding wave function has n zeros, which proves the existence of precisely 
n bound states. The black hole case is considered in Sec. 8, the results apply also for the 
case of charged non-Abelian black holes found in the EYM theory with the SU (2) x f/(l) 
gauge group The second variation of the ADM mass functional for the regular BK 
solutions is explicitly derived in the Appendix. 

2 EYM field equations 

We start from the action of the SU (2) EYM theory 



S-- 

where F^j^ = d^A^ — d^A^ — z[y4^, Aj,] is the matrix valued gauge field tensor, e is the 
gauge coupling constant, = A^r°/2, and r° (a = 1, 2, 3) are the Pauli matrices. 

In the spherically symmetric case, it is convenient to represent the spacetime metric 
in the following form: 

ds" = ll { (1 - —)a^de - r\dd^ + sin2 , (2) 

\ r 1 — 2m I r j 



where m and a depend on t and r. Here the EYM length scale 1^ = vAttG/c is explicitly 
displayed, after such a rescaling, all other quantities in the theory become dimensionless. 
The spherically symmetric SU (2) Yang-Mills field can be parameterized by 

A = WqLi dt + WiLi dr + {p2 L2-(l-pi) Lg} di3+{{l-pi) L2+P2 4}sin^9 d^, (3) 

where Wq, Wi, pi, p2 are functions of t and r, Li = n"r"/2, L2 = d^Li, sini^Lg = d^Li, 
and = (sin-t? cos (/?, sin^? sin 99, cosi?). The gauge transformation 

A^ UAU-^ + iUdU-\ with U = exp{in{t,r)Li), (4) 

preserves the form of the field (||), altering the functions Wq, Wi, pi, p2 as 

Wo^Wo + dtn, Wi^Wi + drn, p±=pl±^p2^exp(±^^])p±. (5) 
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Consider also parity transformation: -(9— s>7r — ^9, y9— + vr. The action of this 
transformation on the fields @, @ is equivalent to the following replacement: 

m^m, a ^ a, pi pi, p2 -p2, Wo -> -Wq, Wi -Wi. (6) 

It is convenient to introduce the complex variable f = Pi + ip2 and its covariant 
derivative D^f = {d^ — iW^)f, as well as the (1 + l)-dimensional field strength W^u = 
d^Wu — duW^, (/i, z/ = 0, 1). The full system of the EYM equations then reads 

d,{r^aW>''') - 2a Im {fD'^ff = 0, (7) 
D,aD^/-^(|/r-l)/ = 0, (8) 

= -^w,.w^^ + ^1 + A^i A/r + ^(i/r - 1)', (9) 

dm = 2N Re DofiDjy, (10) 
9A^^ = -(T^\Dof\' + \DJ\'), (11) 



where = 1 — 2m/r, asterisk denotes complex conjugation. 

In the static case, with Wq = Wi = p2 = and f = Re f = pi, these equations 
are known to possess regular, asymptotically flat solutions discovered by Bartnik and 
McKinnon (BK) |1|], as well as the EYM black hole solutions 0. The regular solutions 
form a discrete family labeled by an integer n = 1, 2, . . .. Near the origin and infinity, 
these solutions, (/n('")5 '"^n(^)c"n('")); have the following common behaviour: 



/ = 1 -6r^ + 0(r^), m = 26V^ + 0(r^), — = 1 + 46V + O(r^), as r 0, (12) 
and 

/ = (-l)"(l-- + 0(l)), m = M + 0(l), a= 1 + 0(4) as r^oo, (^3) 

where < b,(7Q < 1, M and a are numerically known constants depending on n. In 
the domain < r < oo, m{r) and a{r) are monotone increasing functions, while /(r) 
oscillates n times around zero value always staying within the stripe — 1 < /(r) < 1. 

The EYM black hole solutions are distinguished by the two parameters, n and r/^, 
where n is the number of nodes of the function /(r) in the region r > r/j, and r/j G (0, oo) 
is the event horizon size. The boundary conditions at the horizon r = are 

P (T 

/ = A + -x + 0(x2), N = Fx + 0{x^), — = l + ^x + 0(x2), (14) 

F Oh rlF^ 

where 

^=^, P = Wl'l). F = l-Ul-lflrl (15) 

'h 

and < depend on {n,rh)- In the domain < r < oo the behaviour of these 

solutions is qualitatively similar to that for the regular case with the boundary conditions 
at infinity defined by (0). 
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In what follows, it will be convenient to introduce the tortoise coordinate p such 
that dp/dr = 1/aN. In the regular case p runs over the half-line [0, oo): 

p= hO(r^) as r ^ 0, p = r + 2M\nr + 0(1) as r -> cx). (16) 

For the black hole solutions the domain of p is the whole line (— oo, +oo) such that near 
the horizon 

AT = 1 - ^ = Fexp {(^hF^^ + O (^exp (2(rhF^^^ as p ^ -oo. (17) 

3 EYM sphalerons and vacuum to vacuum paths 

As was shown in Ref.0], the odd-ra BK solutions may be treated as sphalerons. This 
interpretation is based on the fact that, through any such a solution, one may find a one- 
parameter family of static field configurations which interpolates between topologically 
distinct EYM vacua. If A denotes the parameter of the family, the corresponding gauge 
field can be chosen as follows 0] 

A^dx'' = ^(1 - Ur))UdU~\ U = expiiXn'^T''), (18) 
where A ranges from zero to vr. The metric functions m and a in (Q) are given by 

(j(r) = exp{—2sin'^ \ [ fn~}^ 
Jr r 

It is worth noting that, for any value of A, the fields of this family satisfy the two 
Einstein equations Gq = SttGTq and Gj; = SnGT^; when A = 7r/2 these equations reduce 
to Eqs.(D, ([n]) |[. One can check that, when A = 7r/2, the fields defined by (|18]), (|19D 
coincide with the n-th BK solution field described above. When A approaches zero or vr 
values, the spacetime metric becomes fiat and the gauge field (|T8|) vanishes. Thus, the 
family of fields (0), (|l^) describes a loop (in the EYM configuration function space) 
which interpolates between the trivial vacuum and the n-th BK solution when A runs 
from to 7r/2, and returns back to the vacuum as A changes from 7t/2 to vr. 

It is important that this loop can be transformed to a path interpolating between 
neighboring EYM vacua with different winding numbers of the gauge field. To see this, 
one has to pass to such a gauge, in which the gauge field potential decays faster then 
1/r as r ^ oo. The gauge transformation 

A^^U{A^ + -d^)U-\ with f/ = exp(z^(/„-l)nV'^) (20) 

yields 

A^dx^ = Y^{1- /n(r))f/+ dU-' + ^(1 + /„(r))t/_ dUz\ 

f/± = exp(zA(/„±l)nV72), (21) 
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whereas the gauge invariant metric functions ([T9|) remain unchanged. One can see that, 
when A runs from zero to vr, the field (0) interpolates between two different pure gauges. 
Notice that, in the new gauge, one has Aa = 0{r) when r — and Aa = 0(l/r^) as 
r — *■ cxo, with Aa being the field component with respect to an orthonormal frame. 

Let A depend adiabatically on time, in such a way that X{t = — oo) = and 
A(t = oo) = vr. Notice that the structure of the field (|2T|) ensures the temporal gauge 
condition Aq = even for time-dependent A. Recall the Chern-Simons (CS) current and 
its divergence 



rtr- 



2ie 



=A^{VaAf3 - —A^Ap), 
9 'J 



V^iT" = — ^^rF^,F'^^ (22) 



where Va is the covariant derivative, and F'^'^ is dual tensor. Integration leads to the 
following relation: 



167r2 



dt / £xy/^trF^^F'"' = / d^x^^K 



+ 



dtl Kdt. 

I J 



(23) 



In the temporal gauge used, the field decays faster then 1/r as r — > oo, such that the 
surface integral entering the right hand side vanishes, while the gauge invariant left hand 
side is H 



dtXsin^Xr dr Wl-l) 



47r^ 



fn 



(24) 



Noting that the gauge potential ( |2TD vanishes at A 
of the CS number along the path (^T]) : 



/„(oo)=-l 

/n(0) = l 

0, one obtains the following value 



sin\cosX){-fl 



CS 



-gK^d^x 



TC 



-(A — sinXcosX), 



(25) 



which changes from zero to one as A runs from to vr. 

Thus the fields (H), ([T9|), (|2TD interpolate between two EYM vacua with different 
winding numbers of the gauge field, and coincide with the (gauge transformed) field of 
the n-th BK solution when X = 7i/2. This allows one to treat odd-n BK solutions as 
sphalerons, "lying" on the top of the potential barrier separating distinct topological 
vacuum sectors of the theory. The profile of this barrier may be obtained from Eq.(p!9D 
due to the fact that the metric function m(r) obeys the initial value constraint, which 
allows one to define the ADM energy [|] [g , 

E[X,f{r)] = lim m(r) = 



sin^X 



if 
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sm X- 



-)exp{—2sin'^ X / /'^ — )dr, 
Jr r 



2r2 



(26) 



such that the sphaleron energy is given by E[tt/2, /n,(r)]. 



4 Perturbations of the equilibrium solutions 



The functional E given by ( p6D is quite useful for a qualitative understanding of the 
behaviour of the potential barrier surface in the vicinity of sphalerons . The sphaleron 
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configurations, A = 7r/2,/(r) = fn{r), correspond to critical points of E[X,f{r)]. One 
can easily see that ^'[A, f{r)] reaches its maximal value at A = tt/2, which shows that each 
sphaleron possesses an unstable mode this kind of instability can be naturally called 
sphaleron instability. In addition, the value of E[7i/2, fn{n)] decreases under certain 
deformations of another type: fn ^ fn + |13[ which will be referred to as gravitational 
instabilities. 

To put the this into a more rigorous form, we return to the general system of the 
EYM field equations ([7|)-([TT|) and consider small spherically-symmetric perturbations of 
a given (regular or black hole) equilibrium solution: 

m^m + 6m, a ^ a + 6a, f^f + 6f, Wq = Wi = ^ SWq, 5Wi. (27) 

In view of (|^), these perturbations can be classified in accordance with their parity. 
Observing that the background equilibrium solutions are invariant under parity, one 
concludes that even-parity perturbations, {6m,6(T,6pi), and odd-parity perturbations, 
{6Wq,6Wi,6p2), must be independent, which is confirmed by the straightforward lin- 
earization of Eqs.(|7D- (p!TD . Notice that the infinitesimal gauge transformation (^ does 
not alter the even-parity perturbations, while the odd-parity ones change as 

6Wo 6Wo + dtn, 6Wi 6Wi + dr^l, 6p2 Sp2 + nf. (28) 

Perturbations in the even-parity sectors have been studied in In this case, 
linearization of Eqs.(^), (p!OD, (p!l|) reveals that perturbations of the metric functions, 
6m and 6a, can be expressed entirely in terms of 6pi, the latter satisfies the following 
Schrodinger-type equation: 

.2 




= uj'(j), (29) 

where 6pi = exp{—iu!t)(j){p) with p being the tortoise coordinate introduced above 
(throughout this paper we use the following notations: y' = dy/dp, and y[. = dy/dr). 
For the n-th background regular or black hole solution, numerical analysis has shown 
the existence of n bound states in (p9|). 

The even parity perturbations correspond to the gravitational perturbations defined 
above. It is instructive to see how the perturbation equation (|29|) can be derived directly 
from the potential barrier functional E |]13[. Let us introduce a harrier height functional 



e[/(r)] = E[A = vr/2,/(r)], (30) 

and consider its expansion near a critical point 

e[f^{r)+<P{r)]=e[Ur)]+6h + ..., (31) 

where the first order term vanishes and dots denote higher order terms. The second 
order term reads 

S'e = 0(-— + V)<f>dp, (32) 
The effective potential V in this expression is explicitly derived in the Appendix, and 



turns out to be the same as that in the Schrodinger equation (p9|). This means that any 
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solution of the perturbation problem (p9|) with u"^ < gives rise to the negative second 
variation of the barrier height functional. 

Let us pass now to the analysis of the odd-parity modes. The odd-parity per- 
turbation sector contains the sphaleron instabilities of the equilibrium solutions. For 
the odd-parity perturbations, the metric remains unperturbed, such that the perturba- 
tion equations can easily be obtained via expanding the Yang-Mills equations (^, (H) 
with respect to SWq, SWi, and Sp2 alone. We use the gauge freedom ( pS]) to impose 
the condition 6Wo = 0, and specify the time dependence as SWi = exp{—iut)a{r), 
Sp2 = exp{—iut)C,{r). The resulting equations read 

f-;f ATaA + _ 1) V + + d^^j^^^ ^ u^^^ ^^^^ 

\ ar ar J ar aJM 

2 

2N^ {fa + fX - fC) = (34) 

where /, a and refer to the background equilibrium solution. There exists also an 
additional equation due to the Gauss constraint (Eq.(^ with z/ = 0): 

-fo=0. (35) 

To simplify these equations we introduce the new variable x = {r'^/2a)a and pass to 
the tortoise coordinate p. Multiplying (|33D by / and using the background Yang-Mills 
equation 

f" = ^fif - 1), (36) 
one can represent Eqs.(p3|-^) in the following form: 

(/Vx + /'e-/e')' = ^Ve, (37) 

fVx + n - fi' = ^'X, (38) 
u^ix' - fO = 0, (39) 
where 7^ = 2Na'^/r'^. One can immediately see that, as long as 7^ 0, only two of these 




equations are independent, say, (|38D and ( pQ]) (one can equally say that Eqs.(|37D, (|38 
are independent, and the Gauss constraint in the case is the differential consequence of 
this equation of motion). When u = 0, Eq.(^) disappears, and there remains only one 
essential equations (|38D which can be solved by any static pure gauge obtained from (p8l): 

X = %, ^ = fn. (40) 



Let u be non- vanishing. Then one can omit the factor lu in ( p9D and specify the following 
independent equations which are clearly equivalent to (^)- (|39|) : 



K' - fi = UW - uj')x, (41) 

X' = /e (42) 
When 0; = 0, Eq.(H3) does not appear in the initial system, however we will retain this 



equation at zero energy too. In this case, this equation (it is sometimes called "strong 
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Gauss constraint") plays the role of the gauge fixing condition. Imposing this condition 

59]) at zero energy, allowing 



violates the remaining gauge invariance of equations (|37|) - 
only for those pure gauge solutions ( |iO| ) whose parameter Q satisfies 




(43) 



5 Transformation to Schrodinger problem 

Our aim now is to reduce Eqs.(|l|), (H) to a one-channel Schrodinger problem. The 
most natural way to do this is to exclude one of the functions, either ^ or x, from the 
equations. In this case however, one obtains either singular or energy-dependent effective 
potential in the resulting Schrodinger equation. For instance, eliminating ^ one finds 



7'(f + l) + 2hV U = ^V, (44) 




where <p = x/ f: the potential in this equation is ill-defined due to nodes of /. Another 
possibility is to try to look for a certain linear combination of x and ^ which might satisfy 
a Schrodinger-type equation with regular, energy-independent potential. Surprisingly, 
this way turns out to be successful. The key observation here is that the combination 

^ = fZx + ^, (45) 

where Z is a solution of the following auxiliary equation: 

Z' = fZ' - 7^ (46) 

satisfies by virtue of (^), (^) and ( PB[ ) the Schrodinger-type equation 

- ^" + U%1) = uj'^i), (47) 

with the energy-independent effective potential 

U = a''N ^^^~^ +2[fz)'. (48) 

The inverse transformation from to the variables x ^i-nd C, reads 

X = ^ {fZij + f'i; - fi;') , e = ^ - fZx- (49) 

In the new potential derived the nodes of / are no more dangerous, such that the potential 
is regular provided that Z is regular. Furthermore, the investigation of the behaviour of 
Z is simplified considerably due the existence of another remarkable relation allowing us 
to represent solutions of Eq. ( ]15| ) in terms of solutions of second-order linear differential 
equation (|43|): 

Z = -QA , where A = ^, (50) 

C+ / PA^dp 
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with C being an integration constant. For each particular fi(p), this correspondence 
provides us with a family of solutions to Eq.(|^) whose members are distinguished by 
values of C in (pOl). However, not all of these solutions are well-behaved. Our aim is to 
specify and C in (^) such that the resulting Z{p) gives rise to well-defined potential 
(||) and transformations (|9|). 

Before passing to the further analysis, one may wonder why it was possible to guess 
the transformation presented. One might suspect the existence of some symmetry in the 
problem, which could naturally explain the appearance of all the relations listed above. 
It turns out that there exists indeed a regular way to proceed in the case, since the 
transformation applied formally resembles the transformation between the dual partners 
in supersymmetric quantum mechanics [T^. First, notice that Eq.(^4D has obvious = 
solutions which are just pure gauge modes (see (^0[)), 



00 = ^'/fl\ (51) 

with Vt satisfying (|13|). Any such solution allows us to factorize the differential operator 
in (H): 



with 



Using (HI) we find 



where 



dp 00 



dp f 



z = -riA, A = -f'^/n'. (55) 

Notice that (^) is equivalent to 

A'/A = fZ, (56) 

and thus Z satisfies the nonlinear differential equation (^61). 

With a standard reduction, the most general solution of the second order linear 
equation (|43|) is 

n = C2A + n (^ci + C2 J^' fA^dp^ , (57) 

where is a special solution, A = and ci, C2 are real constants. This gives 

immediately for Z in (|55|) 

A2 

Z = -QA , (58) 

Cl/C2 + / PA^dp 
Jo 



which is obviously identical to (|50|) . Since the differential equation (|4^) is, for any in 
([57|), identical to 

= cuV, (59) 

one can pass from to 

^ = Q-<p. (60) 
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For C(j 7^ this has the unique inverse 

and by applying Q~ we obtain the "dual" eigenvalue equation 
The differential operator on the left is 



(61) 



(62) 



(63) 



with U being given by (|^), whereas transformations (|60D, ([61| ) are identical to (^ 
(1491). 



6 Structure of the potential 

Now we turn to the investigation of the behaviour of Z. First, consider the case of the 
regular background equilibrium solutions. It is convenient to return for a moment to the 
variable r, then Eq.(^3D takes the form 

2/2 



r'Ql + 2r(l - f;')Q[. 



N 



-n. 



(64) 



This differential equation has regular singular points at r = 0, oo. As is known |T0| , 
the background solutions (p^, (p^) are analytic in some neighborhoods of the origin and 



infinity thus ensuring analyticity of the coefficients in (p^). This guarantees the existence 



of the fundamental system of solutions in the vicinity of each singular point [19|. Using 
one can represent these solutions in the following form: 



n(r) = ci r {1 + i?i(r)} + — {1 + R2{r)} as r - 
and introducing s = 1/r, 
n{s) = rfi {1 + Siis)} + d2-{l + S2is) + A sHn{s)il + Siis))} 



0, 



as s 



(65) 



0, (66) 



where Ci^2 and di^2 are arbitrary constants, -Ri,2('") and Si^2{s) are convergent power series 
vanishing at r = and s = 0, respectively, the quantity A = A{a, M) is specified by the 
asymptotic behaviour of the background solution. 

In the domain < r < cxd Eq. ( 16^) has no singular points and the coefficients of the 
equation are at least continuous functions in this region. This guarantees the existence 
of an extension of the asymptotic solutions (pSf), ( p6| ) to the whole half- line, and this 
extension is at least C^. 

Now we choose a special solution Vt which is regular at infinity. For this we specify 
in (|66D (ii = 1, = 0. Then, near the origin, VL is given by ([651) . One can see that 
the coefficient C2 in this formula does not vanish. Indeed, multiplying Eq.(^) by and 
integrating from p to infinity one obtains 







'or 

.7^ 



+ mdp 



n'. 



+ 



oo 0'2 

{\ + fn^)dp. 



(67) 
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Applying this relation to the particular solution Cl, Cl{oo) = Cl'{oo) = 0, one arrives at 



m=-2^' r{%+fn')dp<o, (6^ 



so that, since Cl vanishes at infinity, it can not vanish also at the origin containing 
therefore the growing branch, ^2 ~ as r — 0. From here we conclude that on the 
domain < r < cxd there exists a smooth (at least C^) solution fl{r) with the following 
behaviour near r = 0, oo: 



fi= ^^ + 0(1) asr ^0, fi = 4 + asr->oo. (69) 

Next, we substitute this solution in Eq.(0), where we require the constant C to 
be positive - to avoid vanishing of the denominator. Then, for each C > 0, the result- 
ing solution Z{r) is smooth and differentiable in the interval (0, oo) with the following 
behaviour near the origin and infinity: 

Z = ao (^^ + . . . + p{Cy + 0{r^)J as r ^ 0, 

2 q(C) ^Mr), 

Z = -- + ... + ^ + 0(-;^) as r^oo. (70) 

Here dots denote terms specified entirely by the background solutions (|12]), ([131) , whereas 
the quantities p{C) and q{C) depend on the constant C as well. Using (|16D we can return 
back to the tortoise coordinate p thus obtaining 

Z=- + 0(p) asp-^0, and Z = -4 + 0(^^) as p ^ oo. (71) 
p p^ p-^ 

We therefore arrive at the one-parameter family of solutions Z{p), whose members are 
distinguished by values of the constant C > in (0), each solution satisfies boundary 
condition ([711). Some of the typical solutions are presented in Fig.l. All these solutions 
have only one zero in the interval (0, oo) - this is because, in view of Eq. (^6[) , the 
derivative Z' is always negative when Z vanishes. One should stress that such a behaviour 
of Z is common for all background regular equilibrium solutions. 

Now, having obtained the function Z, we can determine the effective potential 
U{p) in the Schrodinger equation {^7\). Notice that the first term in ( ^8[ ) contains an 



unbounded piece: cr^A^(3/^ — l)/r^ when p —>■ 0. However, in account of 

condition (^), this divergency is exactly canceled by the second term in (^), such that 
the resulting potential turns out to be continuous in the whole semi-axis, p e [0, oo), 
with the following behaviour near the origin and infinity: 

U{p) = 66(66 - l)al + 0{p) as p ^ 0, U{p) = -| + 0(^^) as p oo, (72) 

p"^ p'^ 

where b comes from Eq . ([T2[) . We therefore find a whole family of regular potentials. 
Each member of this family, f/(p), is specified by the corresponding solutions Z{p) and 
satisfies boundary conditions (ff2[) . A typical example of such a potential is depicted in 
Fig.2. 
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log(p) 

Figure 1: Behavior of the function Z{p) for n = 1 and n = 2 background BK solutions. 
Instead of Z, it is convenient to display Q such that Z = aQ/r. For each n, a pair of 
solutions is shown corresponding to the two different choices of the integration constant 
p{C) in Eq.(70) (the higher maximums correspond to p{C) — 0). 




3 



log{p) 



Figure 2: Odd-parity and even-parity potentials for the n = 1 background BK solution. 
The odd-parity potential is given by Eq.(48), the corresponding Z is shown in Fig.l for 
p(C) = 0. The even-parity potential is defined by Eq.(29). 
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Thus the potential U in the Schrodinger equation (^71) indeed can be chosen to be 
regular everywhere in the half-line [0, oo). However, U is not uniquely defined. Our aim 
now is to show that all potentials from the family specified above are equivalent as far 
as the bound state problem is concerned. To this end we study the properties of the 
transformations (^Sj), (^9]). 

Consider a solution of the initial system of equations (0), (|4^) corresponding to a 
bound state with cu^ < 0. Then, near the origin and infinity, 

^~p, x~P^ asp^O, ^~±-x~ exp{-\Lj\p) as p ^ oo. (73) 

UJ 

Applying transformation (^) with some particular Z{p) and taking (^) into account 
one obtains the solution of Eq.(^7D with the following behaviour: 



'?/'~pasp— >0, ip exp{—\ij\p) as p ^ oo, (74) 
which ensures normalizability of ip 

(iPliP) = ^2 ^p < oo; (75) 
Jo 

so that ip corresponds to the bound state solution of the Schrodinger equation (|^) with 
the same tu^ < 0. On the other hand, starting from (^) and applying the inverse 
transformation (^91), one arrives at (^). It is worth noting that these transformations 
uniquely fix the value ip{0) = ensuring essential self-adjointness of the Schrodinger 
operator H = —d?/dp^ + f/(p) (the Schrodinger equation by itself does not specify 'ipifi) 
because the potential is globally regular in the case). 

Thus, for each particular [/, transformations (^5]), ( ^9]) establish one-to-one corre- 



spondence between bound states with non-zero energy for the initial system of equations 
(|41|), ( ^2]) and those for the Schrodinger equation (|i7|). Obviously, this means that the 
different choices of the potential U are equivalent. Precisely, the nature of this equiv- 
alency is explained by the following considerations. Let Zi and Z2 be two different 
solutions of Eq.(^) (relating of course to the same background field configuration), and 
Ui, U2 are the corresponding potentials. Then, given a solution of ( ^7[ ) for the potential 
Ui and some a;^ ^ 0, ipi, one can construct a solution for the potential U2 and the same 
uj^ due to the following rule: 

^2 = ^1 + ^/(^2 - ^1) (f ^1^1 + /Vi - M) • (76) 

Suppose that V'l is a normalizable solution corresponding to a bound state of \J\, then 
it satisfies ([7^) . From ([70|) , ([16|) one concludes that 



Zi — Z2 ~ p^ as p — s> 0, Zi — Z2 ~ — r as p — » 00, (77) 

P 

which ensures that ?/'2 also satisfies ([71[), that is, it relates to a bound state of f/2 with the 
same eigenvalue up' < 0. Thus the bound state problem for Eqs.(^), (|42D is equivalent 
to that for the Schrodinger equation ( ^7D with any of the effective potentials from the 
family specified above. 
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7 Number of bound states 



Our purpose now is to show that Eq.(^7D admits a zero energy bound state. The simplest 
way to see this is to start from Eqs.(^T]), (^2p which are known to possess the pure 
gauge solutions with zero energy. These solutions are given by Eq.(|y) with the gauge 
parameter Q satisfying Eq.(^). Then, substituting (^Dj) in one obtains a zero energy 
solution for Eq. (^7]) , ipo = f(j), where (j) = Q + ZQ' /'j'^. Next, observing that the function 
(j) obeys the following equation, (p' = pZcp, one can represent ipo as 



2Po = fexp(^J'j'Zdpy {71 



with some positive po', the validity of this solution can also be checked directly. The 
exponential factor in this expression is manifestly positive when < p < oo and, on 
account of the boundary conditions for Z, vanishes at p = 0, oo, such that 

V'o ~ P as p — > 0, ~ ^ as p — > oo, (79) 

P 



which ensures normalizability of ipQ. Thus ipo, being a normalizable solution of (|47|) and 
satisfying the right boundary condition at the origin, ipo{0) = 0, relates in fact to zero 
energy bound state. One should stress that this particular bound state occurs only for 
the Schrodinger equation { ^7\ j and not for the initial system of equations (|1), (|2). 



In fact, the equivalence between the initial system and the Schrodinger equation holds 
only for non-zero energy. When uj"^ = 0, transformation (^) becomes ill-defined, and 
the equivalence is lost. Then, taking into account Eq.(|69D one can see that the pure 
gauge solutions (^Op by themselves can be neither normalizable nor even just regular 
both at the origin and infinity. However, their certain linear combination, i/jq, can be 
normalizable. The idea is to use ipo to establish the existence of negative energy bound 
states for the Schrodinger equation, and those are certainly equivalent to bound states 
appearing in the initial system. 

To this end we observe that, owing to the oscillating factor / in (^), ipo has 
exactly n nodes, which clearly proves that the Schrodinger equation admits n bound 
state solutions with negative energies. We therefore conclude that for the n-th regular BK 
solution there exist precisely n negative modes in the odd-parity spherical perturbation 
sector. 



8 Black hole case 

For black holes, the proof is essentially the same as that in the regular case. We start 
again form Eq. (^3D and want to specify a particular solution for this equation, f^(p), 
such that the corresponding Z{p) is well behaved. The variable p now runs over the line 
— oo < p < oo, such that p = — oo (horizon) and p = +oo are regular singular points of 
the differential equation. Instead of p, it is convenient to use the independent variable x 
defined by ([TSD , then Eq.(|6^ takes the form 

x{l + xf fiL + 2x(l + x) (l - {fS /rl) fi; = 2f^ n. (80) 
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In this equation, the coefficients are analytic functions in the vicinity of x = 0, because 
the background solutions (|I^) are analytic in some neighborhood of the horizon [1^. This 



ensures the existence of the fundamental system of solutions in the vicinity of x = 0. We 
choose the particular solution Q whose regular integral vanishes at the horizon: 

fi(x) = x(l + X(x)), as x-^0, (81) 

where X{x) is given by a power series convergent in the vicinity of x = 0, X{0) = 0. 
Since the coefficients in ( PUD are at least continuous in the domain x > 0, the solution 
([81| ) admits at a smooth extension into this region. When x tends to infinity, the solution 
is given by (|66|) where, repeating the arguments used in Eqs. (|67|) , (168|) , the coefficient 
before the growing part can be shown to be non-vanishing. From here one concludes 
that on the domain < x < oo there exists a smooth solution of (|80| ) with the following 
behaviour near x = 0, oo: 

n = X + O(x^) as X — > 0, and = const x + 0(1) as x — > oo. (82) 

Next, using ([TB|), (|T^) we express this solution in terms of p and substitute it into Eq.(^), 
where, in order to ensure the non-vanishing of the denominator, we choose the constant 
C such that 

C<- fA'dp; (83) 

JO 

the integral here exists since Eq.(^2|) ensures that A ~ 1/p^ when p oo. For each 
certain C we obtain from (^) a solution Z{p) which is continuous and differentiable in 
(— oo, +oo) with the following behaviour near the horizon and infinity: 

^ 0{\) as p-^-oo, and Z -> --^ 0(^^^) as p oo. (84) 



fhP ' >^ 

The corresponding effective potential U{p) is smooth and bounded. The equivalency of 
the different particular U can be checked with the use of (|76D, (^). Then, repeating the 
analysis performed above for the regular case, we arrive at the zero energy solution ipQ 
given by ([78|). In view of (0), this function satisfies the following boundary conditions: 

ijjo ^ - as p —>■ — oo, V'o ~ as p —>■ oo (85) 

p p2 

and is therefore normalizable on the whole line (— oo, +oo). Thus in the black hole case 
there exists a zero energy bound state too. The function ipQ possesses n nodes, which 
shows the existence of n negative energy bound states. From here we conclude that for 
each (n, r/j) EYM black hole solution, Vh G (0, oo), there exist precisely n negative modes 
in the odd-parity spherical perturbation sector. For the n = 1 EYM black holes, the 
negative modes energy eigenvalues for the both parity types are presented numerically 
in Tab.l. 



Tab.l. The bound state energies for the n = 1 EYM black holes. 
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rh 


uj'^ (even-parity) 


cj^ (odd-parity) 


(regular case) 


-0.05249 


-0.06190 


0.1 


-0.05026 


-0.06182 


0.5 


-0.04021 


-0.06000 


1 


-0.02683 


-0.04926 


5 


-0.00212 


-0.00298 


10 


-0.00054 


-0.00075 



Finally, we consider the case of the charged SU{2) x f/(l) EYM black holes ||T^ 
which possess the additional U{1) dyon type field. In this case, the only new thing 
appears - an extra term g^/2r^ in the right hand side of Eq.(^, where = q"^ + q'^, 
qe and q^ being the dyon electric and magnetic charges; the odd-parity perturbation 
equations remain intact. The background solutions are qualitatively the same as in the 
g = case, but the presence of the charge entails now that Vh > q. When rh > q, all 
of the analysis performed above for the q = case remains valid. The only difference 
arises in the extreme limit, Vh = q, when N has zero of second order at the horizon, and 
Pi'^^h) = 1- After all, the change of the boundary condition at the horizon results in the 
following change in the asymptotic behaviour for ipQ given by (|75[): 

■^0 ~ ^ as p ^ —oo, V'o ~ as p ^ oo. (86) 
p'^ p^ 

One can understand this as follows: for extreme black holes the geometry in the limit 
p — ^ — oo is known to be conformally flat. The dynamics of the conformally invariant 
Yang-Mills fleld then must be the same as that in the asymptotically flat region, p — > oo 
(remind that the odd-parity modes involve perturbations of the Yang-Mills fleld alone). 
Thus ipQ is normalizable in the case as well, and we therefore conclude that all charged 
SU (2) X f/(l) black holes, including extreme solutions, also possess n odd-parity spherical 
negative modes. 
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Appendix A Variation of the energy functional 

In this Appendix second variation of the ADM mass functional presented in the main 
text is explicitly derived. The analysis is carried out for the regular BK solution, the 
generalization on the black hole case is straightforward. 

Consider the barrier height functional deflned by Eq. (|30D : 

e[f{r)] = / if" + )exp{-2 / f''-)dr. (Al) 

Jo ZT Jr T 

Let /(r) be a smooth function satisfying the following conditions 

/'(r) = 0(r) asr-^0; /'(r) = ©(l/r^) as r ^ oo. (A2) 
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Define for convenience two new functions 

a{r) = exp{-2 r r-}, m{r) = ^ T (f" + ^J^^)adr, (A3) 
Jr r a[r) Jo r-^ 

they satisfy the following boundary conditions 

a(oo) = 1, (t(0) 7^ 0; m(oo) < oo, m(r) = O(r^) as r 0, (A4) 
and also, by definition, the following equations 

<y' = r—a, {ma)' = (/'^ + ^\ , ' )<y- (A5) 
Consider small variation 

/(r) ^ /(r) + ¥.(r), (A6) 

where 

(^(0) = (^(oo) = 0. (A7) 
To preserve the boundary conditions, the variation must also satisfy 

(y9'(r) = 0(r) as r — > 0, Lp'{r) = 0{l/r'^) as r — > cxo. (A8) 

Putting ( |X6|) into ([Al|) and expanding the result over ip, one obtains 

e[f + ip]=e[f] + 5e + 5h + ..., (A9) 

where the first variation is 



oo 

I I 



6e = 2 I {f^'a + ' ' - 21{f" + '-±—^)a}dr, (AlO) 



(f o -\- 

and the second variation is 



6h = I {V + ^-^-r^^' - 8/(/V + , 

JO r 



+ (8/^ - 2J)(r + ^A^»dr. (All) 

dots in (|A9|) denote higher order terms, and the following new functions have been 
introduced: 

/(r)= T/V-, J(r)= r^'^^. (A12) 

Jr r Jr r 

The boundary conditions ( |A2| ) , (|A8|) imply that 



/(O) < oo, J(0) < oo, /(oo) = J(oo) = 0. (A13) 
Consider the first variation ( |A10| ). Using (^), represent 6e as follows 

Se = 2 {f'ip'a + ^^^^^ — '-ipa - 2/(ma)'}rfr. (A14) 
JO 

Integrating by parts one has 

66 = {2fa^ - Alma) 1^ + 2 / {-(/V) V + (y^^^, — + 2mal'}dr. (A15) 

JO 
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Note, that the boundary terms in this expression vanish. Finding I' from ( A12|) and 
integrating by parts once more, one finally arrives at 



5e = 2 



: r{-((l - —)^f)' + ^^^\ ^^ a]^dr. (A16) 
One can see that the vanishing of the first variation implies the Yang-Mills equation 

j(p-i) 



r 



cr- 



(A17) 



Assume now that the first variation vanishes and consider the second variation. 
Using (|A5|) one obtains 



5\ 



'^{a^'' + - 8/(/V + ^\ )a+ 

r"^ 



+ 8/^(ma)' - 2J(ma)'}rfr. (A18) 
Integrating the fourth term in the integrand by parts, and using Eq.( |A12| ), one obtains 

8 / l'^{ma)'dr = Sl^ma - 16 / Il'madr = 16 / Imafip'—, (A19) 
Jo JO Jo r 



where the boundary terms vanish. Combining this result with the third term in Eq.( [A18| ), 
one has 

- 8 / /{(I - — )a/V + l^L^a^^dr. (A20) 
Using ( A17 ),( |Al2| ),( |A5|) , represent this expression as follows 



-8 r/((l - —)af^)'dT = - 8/(1 - — )a/V 
Jo r 



r 

+ 8 r /'(I - ^)a/Vrfr = -8 H {1 - ^-^)af"^^'^ 
Jo r Jo r r 



+ 



-4 / (1 )a'<^<^'dr, 

Jo r 



(A21) 



where the boundary terms vanish. The fifth term in (|A18|) yields 

rao oo poo roo 

-2 / J{ma)'dr = -2Jma + 2 / J'madr = -2 / maip'2 
Jo Jo Jo 



dr 
r 



if 

-2ma—ip 
r 



+ 2 (ma — )'Lpdr 
Jo r 



,2m 



r 



a({)')'ipdr. 



The first term in (|A18 



IS 



aLp''^dr 



{oip')'dr, 



10 JO 

where the boundary terms are zero. Consider also the following expression 

2m, 



= /"((I _ fHy^^ydr. 

Jo r 



(A22) 



(A23) 



(A24) 



Adding the equations ( |A21|) -( [7^^ ) and introducing also the tortoise coordinate, p, 

dr , 2m, 

- = a 1 , A25 

dp r 

one finally arrives at 

/•CO rp 

= ^{-— + V)v dp, (A26) 

where 2 

V = ail- — ){2K(1 - ^))' + ^^a}, (A27) 

which agrees with the potential in Eq.(p9|) provided that Eqs. (|A5|) , (|A17D are taken into 
account. 



19 



References 

[1] R. Bartnik, J. McKinnon, Phys. Rev. Lett., 61 (1988) 141. 

[2] M.S. Volkov, D.V. Gal'tsov, Pis'ma Zh. Eksp. Teor. Fiz., 50 (1989) 312 
[JETP Lett., 50 (1990) 345]; Sov. J. NucLPhys. 51 (1990) 747; 
H.P. Kunzle, A.K.M. Masood-ul-Alam, J. Math. Phys. 31 (1990) 928; 
P. Bizon, Phys. Rev. Lett., 64 (1990) 2644. 

[3] N. Straumann, Z.H. Zhou, Phys. Lett., B 237 (1990) 353; 
N. Straumann, Z.H. Zhou, Phys.Lett., B 243 (1991) 53. 
Z.H. Zhou, N. Straumann, Nucl. Phys., B 369 (1991) 180; 
Z.H. Zhou, Helv. Phys. Acta, 65 (1992) 767. 

[4] D.V. Gal'tsov, M.S. Volkov, Phys. Lett., B 273 (1991) 255; 

[5] N.S.Manton, Phys.Rev., D 28 (1983) 2019. 

[6] F.R.Klinkhamer, N.S.Manton, Phys.Rev., D 30 (1984) 2212. 

[7] P.Forgacs, Z.Horath, Phys.Lett., B 138 (1984) 397. 

[8] J.A.Smoller, A.G.Wasserman, S.T.Yau, and J.B.McLeod, 
Comm.Math.Phys. 143 (1991) 115. 

[9] J.A.Smoller, A.G.Wasserman,Gomm.Math.Phys. 151 (1993) 303; 
ibid. 154 (1993) 377. 

[10] P. Breitenlohner, P. Forgacs, D. Maison, Comm. Math. Phys. 163 (1994) 141. 

[11] J.Burzlaff, NucLPhys., B 233 (1984) 262. 

[12] M.S. Volkov, Phys. Lett., B 334 (1994) 40. 

[13] M.S. Volkov, Phys. Lett., B 328 (1994) 89. 

[14] M.S.Volkov, D.V.Gal'tsov, Phys.Lett., B 341 (1995) 279. 

[15] O. Brodbeck, N. Straumann, Phys. Lett., B 324 (1994) 309. 

[16] G.Lavrelashvih, D.Maison., Phys.Lett., 343 (1995) 214. 

[17] M.S.Volkov, O.Brodbeck, G.Lavrelashvih, N.Straumann, Phys.Lett., B 349 (1995) 
438. 

[18] D.V. Gal'tsov, M.S. Volkov, Phys. Lett., B 274 (1992) 173; 

[19] E.T.Whittaker, G.N.Watson, A Course of Modern Analysis, Cambridge University 
Press (1927). 



20 



